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Abstract 

We consider the question of diagonal Riccati stability for a pair of 
real matrices A, B. A necessary and sufficient condition for diagonal 
Riccati stability is derived and applications of this to two distinct cases 
are presented. We also describe some motivations for this question 
arising in the theory of generalised Lotka-Volterra systems. 
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1 Introduction and Preliminaries 

We consider the following problem. Given £ W nxn , determine condi¬ 
tions for the existence of diagonal positive definite matrices P. Q in W ixn 
such that 

A t P + PA + Q + PBQ- 1 B t P -< 0, (1) 

where M 0 denotes that M = M T is negative definite. Throughout 
the paper, M >- 0 (M V 0) denotes that M is positive definite (positive 
semi-definite); M A 0 denotes that M is negative semi-definite. 

When diagonal positive definite solutions P, Q of (Hj) exist, we say that the 
pair A, B is diagonally Riccati stable. 

Our interest in the question stems from the stability of time-delay systems. 
Specifically, when such a pair P, Q exists, the linear time-delay system asso¬ 
ciated with A, B admits a Lyapunov-Krasovskii functional of a particularly 
simple form [Ij. The more general question of when positive definite solu¬ 
tions of m, not necessarily diagonal, exist was highlighted in [TB] and some 
preliminary results on this question were also described in this reference. 
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We now introduce some notation and terminology, and recall some basic 
results that will be needed later in the paper. 

A matrix A in W nxn is Metzler if its off-diagonal elements are nonnegative: 
formally, > 0 for i j. For vectors v,w in M n , v > w is understood 
componentwise and means that Vi > Wi for 1 < i < n. Similarly, v > w 
means v > w, v w and means > Wi for 1 < i < n. 

For a positive integer n, we denote by Sym(n, M) the space of n x n real 
symmetric matrices. For A £ M nxn , A T denotes the tranpose of A. A 
matrix A £ M nxn is Hurwitz if all of its eigenvalues have negative real parts. 
It is classical that A is Hurwitz if and only if there exists some P y 0 with 
A T P + PA -< 0. 

The following result recalls some well-known facts concerning Metzler ma¬ 
trices [14]. 

Proposition 1.1 Let A E M nxn be Metzler. The following are equivalent: 

(i) A is Hurwitz; 

(ii) there exists some vector in M n with Av <C 0; 

(in) there exists some positive definite diagonal matrix D E M nxn with 
A t D + DA -< 0; 

(iv) for every non-zero v E M n , there is some index i with vi{Av)i < 0. 

The next result concerning the Riccati equation is based on the Schur com¬ 
plement and follows from Theorem 7.7.6 in jlj. 

Lemma 1.1 Let A E M nxn , B E M nXn be given. 

Q y 0 in M nxn satisfy (Jf) if and only if: 

A t P + PA + Q PB 
B T P -Q 

Theorem 1.1 115) Let C\, C 2 be non-empty convex subsets of a Euclidean 
Space E with inner product (-,■). Further, assume that C -2 is a cone and 
C\ FI C -2 is empty. Then there exists a non-zero vector v E E such that 

(v,x) > 0 for all x £ Ci, 

(v,x) < 0 for all x £ C 2 . 

In this paper, we apply the above result to convex sets in the Euclidean 
space Sym(m ., M) ofmxm real symmetric matrices equipped with the inner 
product 

(A,B) = trace(AH). 



Then the matrices P y 0, 
) -< 0 . ( 2 ) 
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2 Diagonal Solutions to the Riccati Equation 


In a recent paper E3, it was shown that if A is Metzler, B is nonnegative 
and A 4- B is Hurwitz, then there exist P >- 0 diagonal and Q >- 0 satisfying 
(fTj) (Q need not be diagonal). Later, we shall show that it is in fact possible 
to choose Q to be diagonal in this case. 

We first present necessary and sufficient conditions for the pair A. B to be 
diagonally Riccati stable. This result echoes the condition for the existence 
of diagonal solutions to the Lyapunov equation given in [3]: see also C2H3J 
[2] for more recent work on diagonal solutions of the Lyapunov equation. 


Theorem 2.1 Let A,Bg M” xn be given. The following are equivalent. 

(i) There exist P >- 0, Q >- 0 diagonal satisfying (fTj). 

(ii) For every non-zero positive semi-definite 


H = 


H ii 
#12 


H 12 
H-22 


m 


p2nx2n 


with diag(Hn) > diag(H 2 2 ), the matrix 


AH n + BH 


12 


(3) 


has a negative diagonal entry. 


Proof: (i) => (ii): 

If there exist diagonal matrices P >- 0, Q >- 
from Lemma o that the matrix 

A t P + PA + Q 
B t P 

Let a non-zero H € M 2nx2n as in © be given satisfying diag(ifn) > 
diag(Lf 22 )- As R is negative definite and H is non-zero, it follows that 

trac e(HR) < 0. (5) 

Expanding the left hand side of flSJ) we see that 

trace(LTnA T P + H U PA + H U Q + H 12 B T P + Hj 2 PB - H 22 Q) < 0 (6) 

which implies that 

trace(2PAffn + 2PBHf 2 + Q(H\i — # 22 )) < 0. 

As diag(-ffn) > diag(lL 22 ) and Q >- 0 is diagonal, it follows that trace(Q(-Hn — 
H 22 )) > 0 and thus we can conclude that 

tr&ce(P(AF[ii + BH ^)) < 0. 



0 satisfying m, it is immediate 
-Q ) - 1 °- (4 > 
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As P P 0 is diagonal, this implies that the matrix 

AH n + BHj 2 

has at least one negative diagonal entry. 

(ii) => (i): We shall prove the contrapositive. Suppose that there exist no 
positive definite diagonal matrices P , Q satisfying dTJ). Consider the sets 

Ci := { ( ATP + B rp +Q PP ^j\Q,PP0 diagonal € M nxn | , (7) 

and 

C 2 := {H G M 2nx2n | H = H t P 0}. 

Then it is easy to see that C\ and C 2 are convex cones and by assumption 
we have that 

Ci (1C 2 = 0 . 

It follows from Theorem ll.ll that there exists a non-zero matrix H in Sym(2n, M) 
such that 

trace(PbY) < 0 for all X G C 2 
trace (HX) > 0 for all X G C\. 

The first of the above two inequalities implies that HP 0. The second 
implies that 

trace(7fnA T P + H U PA + H l2 B T P + Hj 2 PB + (H n - H 22 )Q) > 0 (8) 

for all choices of diagonal matrices P P 0, Q P 0. Rearranging, we see that 
for all P P 0, Q P 0 diagonal 

trace(2P(AHn + BHf 2 )) > trace(Q(P 22 — -^li))- (9) 

If there exists a diagonal Q P 0 such that trace(<5(1722 — Bn)) > 0, then for 
a given diagonal P P 0, we could obtain a contradiction to Q by considering 
tQ for t > 0 sufficiently large. This implies that trace(Q(i722 — Hu)) < 0 
for all diagonal Q P 0. It follows from this that di&g(H 22 ) < diag(77n). 
Moreover, for every diagonal P P 0, we must have 

trace(P(Affn + BHj 2 )) > 0. 

Now suppose that AHu + BHj 2 has a negative diagonal entry: say in the ith 
position. Then we could choose P P 0 diagonal such that trace(P(APn + 
BHf 2 )) < 0 by setting the ith entry of P to be large enough and keeping the 
other entries of P small and positive. It follows that under the assumption 
that the pair A, B is not diagonally Riccati stable, that all diagonal entries 
of AHn + BHj 2 are nonnegative. Hence, if AHu + BHf 2 has a negative 
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diagonal entry for every H satisfying the conditions of the theorem, there 
must exist diagonal P >- 0, Q P 0 satisfying ©■ This completes the proof. 

A simple application of Schur complements shows that if the pair A, B is 
diagonally Riccati stable and P, Q satisfy (JT]) , then 

A t P + PA -< 0 
(A + B) t P + P(A + B) P 0. 

In the paper m, it was asked whether or not the converse of this statement 
held. It is possible to use Theorem 12.11 to answer this and show that it is 
not the case. 

Example 2.1 Consider 



It is easy to verify that 



satisfies A T D + DA 0, (A + B) T D + D(A + B) P 0. However, if we 
consider H of the form © with 

3 )’ ffl2= ( ~0 -1 )’ H22 = (l “)• 

then HP 0 and diag(H 22 ) < diag(Hn). However, it can be verified by direct 
computation that AH\\ + BH^ 2 has no negative diagonal entries. Hence, by 
Theorem \2.1[ there do not exist diagonal P P 0, Q P 0 satisfying ©■ 

3 Applications 

We now present a number of applications of the general result Theorem l2.ll 
first considering the case when A is Metzler, and B is nonnegative. 

Theorem 3.1 Let A € M nxn be Metzler and B € R nxn be nonnegative. 
There exist diagonal matrices P P 0, Q P 0 satisfying © if and only if 
A A B is Hurwitz. 

Proof: We know that if such matrices P. Q exist, then there exists a diagonal 
Dp 0 such that (A + B) T D + D(A + B) -< 0 and, in particular, that A + B 
is Hurwitz. 

Conversely, suppose that A + B is Hurwitz. Let a non-zero 

*=(w l) do) 
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in Sym(2n,M.) be given with diag(W) < diag(t/). We wish to show that the 
matrix AU + BV T has a negative diagonal entry. To this end, note that as 
HP 0, it follows that U P 0 and hence that 



for 1 < i,j < n and of course un > 0 for 1 < i < n. Furthermore, as H P 0, 
we also have that 



Writing g for the vector in 1" whose zth component is given by y/ufi, g / 0 
as otherwise H = 0. As A is Metzler and B is nonnegative, we have 

n n 

(-AU T BV )a — ^ ( CLijUji T ^ ( bjjVjj 

3 = 1 i =1 

n n 

— ^ Q j ij9j9i d~ ^ ^ ^ijQjQi 

3 = 1 i =1 

= gi[{A + B)g]i. 

As A + B is Metzler and Hurwitz and g ^ 0, it follows from ProDosition ll.il 
that there exists some index i such that gi[(A + B)g\i < 0 and hence that 
(AU + BV T )u . It now follows from Theorem 12.11 that there exist diagonal 
matrices P >- 0, Q y 0 satisfying CD as claimed. 

Our next result gives conditions under which lower triangular pairs A, B are 
diagonally Riccati stable. 

Theorem 3.2 Let A £ M nXTl ; B £ M nxn be lower triangular matrices. Then 
there exist diagonal matrices P >~ 0, Q y 0 satisfying CD if and only if A is 
Hurwitz and \ba\ < \aa\ for 1 < i < n. 

Proof: If there exist diagonal matrices P y 0, Q >~ 0 satisfying CD then 
A t D+DA -< 0 and hence A is Hurwitz. To show that \bu\ < \au\ for all i. we 
argue by contradiction. Assume that there is some i with \ba\ > \aa\. First 
assume that bn > 0 (the case bn < 0 can be handled identically.) Consider 
the matrix H £ Sym(2n, M) of the form (HOD where 

U = V = W = e ie J. 

So U and V are both the diagonal matrix with all entries zero apart from 
the zth diagonal entry. It is then easy to see that HP 0 and that diag(W) < 
diag(17). The only non-zero diagonal entry in AU + BV T is 

(AU + BV t )u = 


> 
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Q'ii'U'ii biiVa 
(Lii bii 

0 


as bu > 0 and \bu\ > \au\. It follows from Theorem 12.11 that there cannot 
exist diagonal matrices P P 0, Q P 0 satisfying (P) which is a contradiction. 
For the case where bn < 0, choose H with U = W = , V = —. 

For the converse, assume A is Hurwitz and that |6jj| < |ajj| for 1 < i < n. 
Consider HP 0 in Sym(2n, R) given by (fTOli with diag(TU) < diag(U). As 
H is non-zero, let i be the lowest index such that un >0. As U is positive 
semi-definite, it follows that u yi = 0 for all j < i. Moreover, from our 
assumptions on H, we must have 



= 0 


for all j < i. Moreover, we also have \vn\ < \un\. Using these observations 
and the fact that A, B are lower triangular, we see that 

n n 

i^AU H - BV ^ ^ Q'ij'Uji + ^ ^ ^ijVij 

3 = 1 3= l 

= ^2 CLijUji + ^2 bijVij 

j<i j<i 

— Q'ii'U'ii T bnVn . 

The result now follows from the simple observation that anun < 0 and 

An identical argument where we choose the index i to be equal to the largest 
index for which u tl > 0 will show that a corresponding result also hold for 
upper triangular matrices. 

Comment: 

Note that Example 12.11 shows that even when both A and B are lower 
triangular and A and A + B are both Hurwitz, it is possible that the pair 
A, B is not diagonally Riccati stable. 

4 Metzler and Nonnegative Matrices - A Closer 
Look 

In the previous section, we showed that if A is Metzler, B nonnegative and 
A + B Hurwitz, there exist diagonal matrices P P 0, Q P 0 satisfying ©• 
We have also noted that for any such pair P, Q, the matrix P must be a 
simultaneous diagonal solution of the Lyapunov inequalities associated with 
A, A + B. We next show that for any diagonal P satisfying A T P + PA -< 0, 
(.A + B) t P + P(A + B) -< 0, there must exist a diagonal Q such that P, Q 
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satisfy Q. Note that in this case, as B is nonnegative, from ( A + B) T P + 
P(A + B) -< 0 it follows that A T P + PA -< 0. 

This strengthens the result given in m and moreover, the method of proof 
used here is significantly different. Moreover, the method of proof gives us 
a constructive approach to finding P and Q. 


Theorem 4.1 Let A £ M nxn be Metzler and B € M. nxn be nonnegative and 
assume that A + B is Hurwitz. Let P 0 be diagonal and such that 

(A + B) T P + P(A + B) ~< 0. (11) 

Then there exists Q >~ 0, diagonal such that P, Q satisfy (HD- 


Proof: As P >~ 0 is diagonal it is straightforward to verify that ( A + 
B) T P + P(A + B ) is again Metzler and, moreover as it is negative definite, 
it is Hurwitz. From Proposition 11.11 there exist vectors v 0, w 3> 0 such 
that 

({A + B) t P + P{A + B))v = -w. (12) 

Now choose a diagonal Q >- 0 such that 


Qv 


B 1 Pv + -w. 


A simple calculation now shows that 


/ A t P + PA + Q PB 
V B t P -Q 



-w/2 \ 
-w/2 ) ' 


(13) 


(14) 


Noting that 


f A t P + PA + Q PB \ 
V B T P -Q ) 


is Metzler, it follows from (1141) that M is Hurwitz. However as M is symmet¬ 
ric this implies that M -< 0 and hence that P, Q satisfy CD- This completes 
the proof. 


5 Generalised Lotka-Volterra Systems 

In this section, we highlight how the results described above can be applied 
to a generalised class of Lotka-Volterra models occurring in population dy¬ 
namics. Specifically, we consider the following generalized Lotka-Volterra 
model of population dynamics in an n species community: 


Xi(t) = gi(xi(t)) 


n n 

c J+ E a vfj OlHO)+ 

3 = 1 i =1 



(15) 
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for 1 < i < n. 

For systems without time-delay, similar models were considered in [8], while 
more recently a switched system version without time-delay was studied in 
m®- Systems of the form (1151) with time delay and with some special forms 
of functions fi were studied, for example, in mm 

Here Xi(t) is the population density of the i-th species; the functions g, : 
[0, +oo) —» [0, +oo) and fi : [0, +oo) —>• [0, +oo) possess special properties 
described below and allow for a generalisation of the classical Lotka-Volterra 
model; Cj, a^-, bjj are constant coefficients. The coefficients c* characterise 
the intrinsic growth rate of the i-th population; the self-interaction terms 
Pu9i(%i)fi(xi ) with pa < 0 reflect the limited resources available in the 
environment; the terms Pijgi(xi)fj(xj) for j / i describe the influence of 
population j on population i. 

Let M" be the nonnegative cone of M n : 

R£ := {x € M n | x > 0}. 


Also, intR” is the interior of M+. 

We make use of some standard results and concepts from the Lyapunov 
stability theory of Functional Differential Equations [Lj. In particular, for 
a given real number r > 0 , C([—r, 0 ], int R”) denotes the space of con¬ 
tinuous functions (j>(6) : [—r, 0 ] —» intMT with the uniform norm ||^|| T = 
max 0 g [_ Ti o] ||^(0)||, and || • || denotes the Euclidean norm of a vector. We 
assume that initial functions for (fl5l) belong to the space C([—r, 0], intR”). 
For a solution x(t) of (fl5l) . we denote by xt the element in C{[— r, 0], int ) 
given by xt(u) = x(t + u) for u € [—r,0]. For (j) i n C([—r,0],intR”), (f)j 
denotes its jth component. 

We consider functions fi and gi satisfying the following properties, which 
are consistent with the standard assumptions made in m® and elsewhere. 
Moreover, these assumptions are clearly satisfied by the classical Lotka- 
Volterra model [9]. 

(i) gi and fi are continuous on [ 0 , +oo). 

(ii) For each initial condition in intM™, there exists a unique, globally 
defined solution to (|15|) . 

(hi) 9i(0) = fi( 0) = 0 , gi(xi) > 0 , fi(xi ) > 0 for Xi > 0 . 

(iv) fi is a strictly increasing function on [ 0 , +oo), and fi{xf) —> +oo as 
Xi —» +oo. 

( v ) /i + °° fi(C)/gi( C) d( = +oo; 

( vi ) fo 1 /ffi(O d C = +°°- 
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The assumptions on the system imply that int M" is an invariant set for (USD. 
For biological reasons, we will consider (1151) with respect to the state space 
int M”. The following result demonstrates the relevance of diagonal Riccati 
stability to Lotka-Volterra systems and provides a practical motivation for 
the question considered in this paper. 

Throughout the remainder of this section, A and B denote the matrices 

A = 1 , B = 

Theorem 5.1 Consider the system described by (|15l) . Assume that: 

1. there exists an interior equilibrium point x; 

2. there exist diagonal matrices P y 0, Q y 0 satisfying ©• 

Then the equilibrium point x is globally asymptotically stable in int M” for 
any value of the delay r. 

Proof: Denote f(x) = (/i(.xi), ..., f n (x n )) T , c = (ci,... ,c n ) T . 

By assumption, x is an equilibrium point in intR” so it satisfies: 

c+(A + B)f(x) = 0. 

Using this, the system m can be rewritten as follows 


/ n 

Xi{t) = gi(xi(t)) \^2a ij (fj(x j {t)) - fj(xj)) 

\j=i 


+ XJ hjUAxjit - r)) - fj{xj)) j , 


i = 1,..., n. 


(16) 


We show that for the system (USD, there exists a Lyapunov-Krasovskii func¬ 
tional of the form 


V&) = E 


rm 


Pi 


i=1 


fi {0 - fi(Xj 
9i{ 0 


-d( + E h] 

3 = 1 


ifMo( u )) - fj( x j)) 2 du, 

(17) 


where pi and jij are positive coefficients. 

In fact, we know by assumption that there exist diagonal matrices P y 0, 
Q y 0 such that 


M = 


A t P + PA + Q PB 
B T P -Q 


A - 2 pi 


for some /3 > 0. It follows readily that for x, y € M n , 

x T PAx + x T PBy + 1/2 {x T Qx — y T Qy) < —f3(x T x + y T y ). (18) 
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Choosing pi, ...,p n and g\ , ...,p n so that P = diag(pi,... ,p n ), Q = diag(2/ii,.. . ,2p n ), 
we claim that V given by (1171) defines a Lyapunov-Krasovskii functional for 
the system (USD- 

First of all, it follows from the assumed properties of the functions fi and 
gi that 

Xi(MO) -®ll) <V(tp) < X2(\W~x \\ T ), 

where xi(C) and X2(C) are Khan functions [7|, and V{<p) —>• +oo as ip(0) 
tends to the boundary of intM" . 

Differentiating the functional (|17D with respect to the system (1161) . we obtain 

n n 

V{x t ) = ^Pi(fi(xi(t)) - fi{xi))^ j a ij (f j {x j (t )) - fj(xj)) 

i= 1 j= 1 


+ ~ fi ( x i )) bij (fj (Xj (t - T)) - fj(Xj)) 

*=1 3 =1 

n n 

+ -^Vj{fj(xj(t - r)) - fj(Xj)) 2 . 

3 =1 3 = 1 

Writing this more compactly, we have 

K(x t ) = (/(®(t)) - f(x)) r PA(f(x(t)) - f(x)) 

+(f(x(t)) - f{x)) T PB{f(x(t - t)) - /(*)) 

+J(/(®(*)) - f( x )) T Q(f( x (t)) ~ fix)) 

- t)) - f {x)) T Q(f (x(t - t)) - f(x)). 

It follows immediately from (1181) that 

n 

V(x t ) < -/3^2{{fi{xi(t)) - fi{xi)) 2 + ( fi{xi(t - t)) - fi{xi)) 2 ) . 

i =1 

This completes the proof. 

Comment: The previous result shows that if a generalised Lotka-Volterra 
system has an interior equilibrium point and if there exist diagonal solutions 
P^0,Q^0to(H]) then the equilibrium is globally asymptotically stable 
independent of delay. Moreover, the stability is established using a simple 
diagonal Lyapunov-Krasovskii functional. We next describe the implication 
of this result for so-called mutualistic systems in which the matrix A is 
Metzler and B is nonnegative. 
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Corollary 5.1 Consider the system (USD- Assume that A is Metzler, B 
is nonnegative and that c > 0. Furthermore, assume that A + B is Hur- 
witz. Then there exists an equilibrium point x in intK” which is globally 
asymptotically stable for all values of the delay r. 

Proof: As A + B is Metzler and Hurwitz and c 3> 0, it follows that 

— (A + B)~ l c 3> 0. 

It now follows immediately from property (iv) of the functions /i, 1 < i < n, 
that there exists a unique x € int HA satisfying 

f(x) = -(A + B)~ l c , 

which implies that x is an equilibrium of (|15|1 . Moreover, as A + B is Metzler 
and Hurwitz, Theorem 13.11 implies that there exist diagonal matrices P y 0, 
Q y 0 satisfying ©• The result now follows from Theorem EHJ 

To finish off this section, we note that the existence of diagonal P, Q sat¬ 
isfying the Riccati inequality also implies uniform boundedness of the gen¬ 
eralised system (|15j) in the case where there is not an interior equilibrium 
point. 

The system (1151) is uniformly ultimately bounded with the ultimate bound 
R > 0 if, for any M > 0, there exists T = T(M) > 0, such that ||x(t, cf>, to)|| < 
R for all to > 0, \\4>\\ T < M, t>to + T. 

Choose an interior point, not necessarily an equilibrium, x € intM”. Con¬ 
sider the Lyapunov-Krasovskii functional (1171) . 

As in the proof of Theorem 15.11 we obtain that V(<f) —>• +oo as <f>{ 0) tends 
to the boundary of int M” . 

A relatively straightforward calculation shows that for the derivative of the 
functional (fT71) with respect to the system (fT5l) the estimate 

n 

V[xt) < ~/3^2{{fi{xi(t)) ~ fi(xi )) 2 + (fi(xi(t - t)) - fi(xi)) 2 ) 

i= 1 


+y^2\fi(xi(t)) - fi(xi )| (19) 

2—1 

holds, where f3 and i] are positive constants. 

This implies that there exists a number L > 0 such that V(xt) < 0 for 
||x(t)|| > L, x{t) € M", which in turn implies that the system (1151) is uni¬ 
formly ultimately bounded in M" . 
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6 Concluding Remarks 


We have considered the general question of when a pair of matrices A, B 
is diagonally Riccati stable. Building on the result given in dj, we have 
derived a characterisation of diagonally Riccati stable pairs the echoes the 
classical result of Barker, Berman and Plemmons for the Lyapunov equa¬ 
tion. Applications to Metzler and Triangular matrices are presented and the 
relevance of diagonal Riccati stability to generalised Lotka-Volterra systems 
in population dynamics has been highlighted. 
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